Mon. Not. R. Astron. Soc. 000, dini (2004) Printed 2 February 2008 (MN WF^ style file v2.2) 



On the incidence of KHz quasi-periodic oscillations in the 
neutron star system 4U 1636-53 

R. Misra^* and K. Shanthi^f 

^ Inter- University Centre for Astronomy and Astrophysics, Post Bag 4, Ganeshkhind, Pune-411007, India 
^ UGC Academic Staff College, University of Mumbai, Mumbai-400098, India 

Accepted - . Received -; in original form 

ABSTRACT 

Nearly 1.32 x 10® sec of non-zero count rate data were obtained from five years of 
RXTE observations of the atoll source 4U 1636-53. The data was divided into 10309 
segments of 128 sec each. The histogram of the number of segments as a function of 
count rate shows that the system can naturally be classified into four fiux states. For 
each segment an automated search for kHz QPO was undertaken and the histograms 
of the maximum Leahy normalised power Pmax for the four fiux states were created. 
These were fitted by probability distribution functions under the assumption that the 
intrinsic amplitude of the QPO has a normal distribution. The best fit distribution 
functions, showed that the fraction / of the time the system exhibits a kHz QPO, 
decreases from near unity for the lowest flux state {fv 1400 c/s) to zero for the highest 
one (> 2550 c/s) while the average amplitude remains nearly constant. Based on 
the best fit probability distribution functions, a threshold Pthre was defined, such 
that 95% of the segments with Pmax > Pthre correspond to a real QPO signal. For 
segments selected by this criterion, the frequency versus count rate plot reproduced 
the earlier known "parallel track" variation with the tracks coinciding roughly with 
the flux states. It is shown that the gaps between the tracks are not caused by uneven 
sampling, but rather the QPO phenomenon is absent or weak when the system's flux 
level is intermediate between two flux states. 
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1 INTRODUCTION 

Rapid (300-1300 Hz), nearly periodic variability in the X- 
ray light curves of low mass X-ray binary systems (LMXB's) 
have been observed by the Rossi X-ray Timing Explorer 
(RXTE) in over two dozen neutron star bearing LMXB's 
(see van der Klis 2000 for a review). These oscillations, re- 
ferred to as kilohertz QPO (quasi-periodic oscillations), of- 
ten have high quality factors (Q=FWHM/frequency), and 
tend to be seen in pairs, with nearly constant frequency 
separation between the two peaks. The high frequencies of 
kilohertz QPO strongly suggests that they are related to 
phenomena taking place in the inner regions of accretion 
disks surrounding neutron stars, making them potentially 
valuable probes of strong gravity and the behaviour of mat- 
ter in such environments. 

Several theoretical ideas have been proposed to explain 
the phenomenology of kilohertz QPO. In these models, one 
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of the frequencies is generally identified as the Keplerian 
frequency of a characteristic radius of an accretion disk (for 
e.g. the innermost orbit). In the sonic point model (Miller, 
Lamb & Psaltis 1998) the second frequency is identified as 
the beat of the primary QPO with the spin of the neutron 
star. Stella & Victri (1999) have proposed a general rela- 
tivistic precession/apsidal motion model wherein the sec- 
ondary QPO is due to the relativistic apsidal motion of this 
characteris tic orbit. On the other hand, in the two oscilla- 
tor model lOsherovich fc TitarchukH (119991) . the secondary 
frequency is due to the transformation of the primary (Ke- 
plerian) frequency in the rotating frame of the neutron star 
magnetosphere. Recently iLee et al~l (l2004ft have proposed 
that the QPO occurs at that radius of the disk where local 
resonance occurs between the Keplerian and epicyclic fre- 
quencies. These models are based on their predictions for 
the variation of the frequency separation with the kilohertz 
QPO frequency and the variation of the kHz QPO frequency 
with that of other QPO observed in a source. However a 
clear consensus regarding which is the viable model has not 
arisen perhaps because most of these dynamic models gen- 
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erally do not address the direct origin of these oscillations. 
Using the observed lags in different energy bands, Lee, Misra 
& Taam (2001), showed that the QPO is driven by a tem- 
perature variation in a Comptonizing cloud of size « 10 km. 
A unified model for kilohertz QPO, that takes into account 
both the dynamic and radiative aspects of the phenomenon 
is still illusive. Such a model should also address how often 
and under what circumstances does a QPO occur and hence 
the primary motivation of this work is to quantify the inci- 
dence of the kHz QPO phenomenon in the atoll source 4U 
1636-53. It is important to know whether the phenomena is 
a ubiquitous or a transient one, which in turn will provide 
insight into its cause. 

In the next section, we argue that the standard data 
analysis technique is not adequate to address this question 
and motivate the need for a modified technique. In that sec- 
tion we highlight the reasons why the atoll source 4U 1636-53 
was chosen for the analysis and the need to study the inci- 
dence in terms of the long term flux states. In §3 the data 
analysis technique and the theoretical probability distribu- 
tion functions which are fitted to the results of the analysis 
are described. In §4 the results are presented and fitted while 
in §5 the implications of the results are discussed. 



2 DATA ANALYSIS TECHNIQUE 

The standard technique to analyse the kHz QPO phe- 
nomenon is to compute power spectra of segments of the 
light curve and fit Lorentzian profiles to them. This pro- 
vides unambiguous measure of the QPO frequency, width 
and amplitude. Correlation between these quantities and 
the spectral properties like intensity, hardness ratio can be 
tested and quantified. However, it is difficult to use this tech- 
nique to obtain an estimate of the incidence of the QPO 
phenomenon. To apply this technique, one has to impose a 
threshold significance level, say 95% confidence level, that a 
QPO has been detected for a given segment of data. This 
means, by definition, that a certain fraction (5% in this ex- 
ample) of the segments analysed will have positive but false 
detections. These detections may be a significant fraction of 
the real QPO that can be detected in this scheme ( say 10%). 
Increasing the threshold level, will decrease the number of 
QPO detected and biases the result toward high amplitude 
detections. Another related problem with this analysis tech- 
nique is the bias in the measure of the amplitude. To illus- 
trate this, consider a long light curve which has a QPO with 
a constant amplitude A, which is less than the threshold sig- 
nificance level of the analysis. The resultant power in each 
segment due to this QPO will not be a constant but will 
have statistical fluctuations. Hence, occasionally, the power 
will be greater than the threshold and this result can be in- 
terpreted as being due to a rarer but stron ger QPO in the 
light curve. As noted bv lLeahv et al. I il983l) . significance of 
the presence and amplitude of a detected QPO, should take 
into account the number of segments and frequency range 
searched. The results obtained without taking these effects 
into account maybe biased and should be interpreted with 
caution. 

These deficiencies of the standard data analysis tech- 
nique are difficult to overcome and hence to make progress, 
signiflcant trade off in the accuracy and interpretation of 



the data have to be made. The technique applied in this 
work, which is described in more detail below, involves find- 
ing the Leahy normalised power spectrum for each segment 
and identifying the maximum power and the frequency at 
which that occurs. The histogram of the maximum power 
can then be compared with theoretical probability distri- 
butions. The advantage of this method is that there are no 
thresholds for QPO detections imposed and hence the data is 
used maximally. However, there are two related drawbacks. 
First, the peak of the power spectrum is representative of 
the peak of the Lorentzian fit to the QPO and not the am- 
plitude. To obtain the amplitude one needs to know also 
the width or Quality factor of the QPO, which is not ex- 
tracted in this analysis. Second, the theoretical probability 
distribution, in principle, should also depend on this width 
or quality factor. Moreover, these probability distributions 
for Leahy normalised maximum power are computed on the 
assumption that the the count rate is a constant for dif- 
ferent segments, which is not true even when the data is 
selected for a given intensity state. Taking these effects cor- 
rectly into account is complicated and hence in this work we 
make simplistic assumptions when computing the theoreti- 
cal probability distributions and partially justify our choice 
by showing that they fit the observed distributions with a 
reduced or order unity. This means that more correct 
but complicated probability distribution functions, perhaps 
cannot be differentiated by the present data and hence may 
not be warranted. 

The presence and properties of the kHz QPO are ex- 
pected to be correlated to the source's spectral parameters 
(see van der Klis 2000 for a review and references therein). 
Low mass X-ray binaries containing weakly magnetised neu- 
tron stars are divided into two classes based on the Z or 
atoll like shaped tracks they trace out in colour-colour di- 
agrams. The time-scale to trace out the tracks is typically 
short (^ day), although atoll sources can take longer (^ 
week) to cover the entire curve. The variation in the inten- 
sity during such a trace out is typically a factor ^ 2, while 
the changes in the colours (depending on the definition) is 
typically ^ 20%. The kHz QPO properties are known to 
be correlated to the position of the source in the colour- 
colour diagram in these timescales. Over longer time-scales 
the behaviour of the system is more complex. For some of 
the atoll sources, there are secular variation in the intensity 
by « 20% while the response corrected colour-colour track 
remains nearly invariant (Di Salvo, Mendez & van der Klis 
2003). The QPO frequency versus intensity plots show par- 
allel tracks where each track marks the short term corre- 
lation between the two (e.g. Mendez et al. 1999 ; Mendez, 
van der Klis & Ford 2001). This perhaps implies that the 
properties of the QPO are determined by the spectral shape 
(i.e. the position in the colour-colour diagram) and not on 
the long term intensity variations. This is partially sup- 
ported by the plot of QPO frequency versus position in the 
colour-colour plot, where the p arallel tracks in the in tensity 
plot almost overlap each other (iDi Salvo et al. l)2003h . How- 
ever, the scatter of the frequencies with respect to colour- 
colour position is larger than the scatter along each track 
in the intensity plot. Moreover, for the long term analy- 
sis there is no significant c orrelation between the frequency 
and colour-colour position JPi Salvo et al. Il2003h . Evidence 
that the QPO properties depends on the long term inten- 
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sity, comes from the frequency-intensity parallel tracks for 
different sources, where for all of them the low intensity 
tracks have a wider range of observed frequencies compared 
to the high intensity tracks. As described in this work, the 
histogram of the number of detections as a function of in- 
tensity shows peaks which nearly coincide with the parallel 
tracks, which implies that the source has distinguishable in- 
tensity (or flux) states. In summary, it is not clear if all of 
the QPO properties depend only on the spectral shape. The 
long term intensity may also be an important factor deter- 
mining some properties especially the incidence. 

In this work we quantify the incidence of the QPO 
phenomenon as a function of different long term intensity 
states. This source was chosen since a highly coherent kHz 
(sa 800 Hz) QPO ha s been observed in its po wer spectrum 
IZhang et all Jl996D : IWiinands et al. I iIIQQTD . The second 
QPO at « 1200 Hz is weak and in general cannot be de- 
tected by the techniques used in this work. 



3 DATA ANALYSIS AND PROBABILITY 
DISTRIBUTION FUNCTIONS 

Nearly all of the archival RXTE AaXa. for the atoll source 4U 
1636-52, that had the Event mode data (E_125^iS_64M_0_ls) 
with 128/isec time resolution, were analysed. The resulting 
1.32 X 10^ sec of non zero count data were divided into 10309 
segments of 128 sec each. For each segment a power spec- 
trum was obtained at a Nyquist frequency of 4096 Hz, by 
averaging 256 sub-segments of 0.5 sec each. Thus the power 
spectra had 2048 frequency bins of width 2 Hz. The power 
spectrum were Leahy normalised jLeahv et al. lll983l i.e. 

P(/) = 2|a(/)|ViV (1) 

where a(/) is the Fast Fourier Transform and A'^ is the total 
number of counts. For each power spectrum the maximum 
power Pmax within the frequency range 200 to 2000 Hz and 
the corresponding frequency fq were obtained. For each seg- 
ment we define a normalised count rate Rn = R{5/npon) 
where R is the detected count rate and Upon is the number 
of proportional counter units (PCU) on during the obser- 
vation. It should be noted that since in the computation of 
Rm, the time dependent detector response is not taken into 
account, it is at best an approximate representation of the 
real flux from the source. 

In the absence of a signal, the average power at any 
frequency is < P 2 and the distribution of P, Pns{P) is 
the distribution with 2M degre es of freedom, where M = 
256, is the number of sub-seements lLeahv et al. I(ll983t ). The 
probability density of obtaining a Pmax is 

Pns (Pmax , ) — IV f Pnsi^Pmax^Q ns(^^ Pmax^ ^ (2) 

where Nf — 900 is the number of frequency bins sampled 
and 

Qnsi<Pmax)^ [ Pns{P)dP (3) 

Jo 

is the probability that the power in a frequency bin is less 
than Pmax- In the presence of a signal the computation of 
the probability density of Pmax, Ps{Pmax) is more compli- 
cated and model dependent. For a pure sinusoidal signal, the 



averag e power at the oscillation frequency is iLeahv et al. I 

<P>=2[l + 0.19Aryl^] (4) 

where A is the amplitude of the oscillation. It has 
been assumed here that the oscillation frequency is much 
smaller than the Nyquist one. The probability distribution 
Ps (P, N, A) can be approximated to be a distribution of 
2M degrees of freedom but rescaled such that the average 
power is given by eqn (4). The total number of counts A' is 
not a constant even for all segments considered in the same 
flux state, since the number of PCU that are on during a 
observation varies. This variation can be taken into account 
by representing the probability distribution as 

5 

Ps{P,A) = J2gMP,N.,A) (5) 

where i is the number of PCU on, g(i) is the fraction of the 
time that i PCU were on and Ni is the average number of 
counts when i PCU were on i.e Ni = RN{npon/5) x 128 sec. 
The above expression assumes that the amplitude of the os- 
cillation is a constant which naturally may not be the case. 
Moreover, in a deflned flux state, the value of Pjv varies over 
a range which will affect the predicted average power (eqn 
4). Further, the observed QPO is not a pure sinusoidal signal 
and has a finite width. Taking these effects into account is 
complicated and hence for simplicity we convolute the prob- 
ability distribution ps (P, A) with a Gaussian of width a to 
give a averaged probability distribution < ps{P,A) >. The 
probability density for obtaining Pmax under these assump- 
tion becomes 

Ps (Pmax , ^) — P ns(,Pmax , N j-' l)C^s(<C Pmax^ ~\~ 

<Ps{Pmax,A) >Q„s{< Pmax ff'\6) 

where 

Qs{<Pmax)= <Ps{P,A)>dP (7) 

Jo 

is the probability that the power in the oscillation frequency 
bin is less than Pmax . The first term in eqn (6) is the proba- 
bility density that Pmax is obtained by chance at any other 
frequency, while the second term is for the probability of 
Pmax to occur at the oscillation frequency. At a given flux 
state, in general, a signal may exist only for a fraction / of 
the time. Hence the total probability distribution is 

P = fPs{Pmax,a,A) + (1 - f)Pns{Pmax,Nf) (8) 

The probability distribution P has three unknown parame- 
ters /, A and a and depends on the known normalised flux 
Rn and PCU on fractions gi. 

The technique used in this analysis is sensitive only 
to narrow band, peaked signals. Any broad band compo- 
nent will be spread over all (or a significant fraction of) the 
900 frequency bins used. Thus the contribution of a equally 
strong broad band signal at a given frequency bin would be 
a large factor (« 100) less than that of a narrow band signal 
which is spread over fewer (« 10) frequency bins. Thus, to 
effect the analysis a broad band signal would need to have 
at least 10 — 100 time more power than a typical QPO de- 
tected (~ 5%). For an average rate of 2000 counts/sec, the 
Poissonian noise at 500 Hz ( 300 Hz) is 100% (30%) and 



4 R. Misra and K. Shanthi 



hence the Poissonian noise is expected to dominate over any 
underlying broad component which justifies the use of 
distribution for pns{P)- However, it should be noted that 
the effect of any such hig h frequency broad band noise (e.g. 
Ivan Straaten et al. 112000) has not been taken into account 
in this analysis. 



4 RESULTS 

The histogram of the normalised count rate Rn for all 
the data is shown in Figure 1. i?jv ranged from ~ 500 
counts/sec to k. 4000 counts/sec, not including the burst pe- 
riods. We have classified the data into four flux states which 
are marked in Figure 1. Data for Rm < 1100 counts/sec 
have not been shown in the figure, since the low count rate 
and the relatively small number of detections does not al- 
low for any significant results to be obtained as discussed 
below. Figure 2 shows the histograms of Pmax for the four 
states. The error for each bin was taken to be the square root 
of the number of detections at that bin. Fitting was done 
over a continuous range of Pmax where there were detec- 
tions. The dotted lines represent the corresponding curves 
if there were no signals in the data PnaiPmax)- Note that 
there are no free parameters in PnaiPmax) and the good fit 
to the high state (IV) histogram (reduced ~ 0.6) implies 
that there are no strong oscillatory signals in this state and 
the fraction of time a QPO is present during this state is 
nearly zero. The upper limit on this fraction / < 0.06 for 
A > 0.025 can be inferred. The shape of the histograms for 
the other three states are clearly different, indicating that 
there are oscillatory signals at least for some fraction of the 
time during these states. These histograms were fitted by 
the probability distribution P for different values of fraction 
/, and amplitude A while a = 0.15 was held constant. The 
best fit curves (solid lines), the best fit parameter values and 
the number of degrees of freedom are shown in in Figure 2. 
Varying a does not change the best fit parameters signifi- 
cantly. For example, for a = 0.25 the best fit values of / 
are 0.97, 0.36 and 0.15 for states I, II and III respectively 
while the corresponding values for a = 0.05 are 0.8, 0.30 
and 0.11. For Rn < 1100 counts/sec, the best fit parame- 
ters obtained were / = l.Otp 4 and A = .0475 ±0.005 which 
are consistent with the interpretation that f ^ 1 for low 
count rates. The theoretical probability distributions ade- 
quately represent the data with reduced ~ 1 for all the 
fits. Hence more sophisticated probability functions are per- 
haps not warranted and may lead to over modelling. 

Based on the best fit probability distribution. We can 
now define a threshold power Pthre such that for 95% of the 
segments with P > Pthre, the peak power corresponds to a 
real signal. The corresponding thresholds for the states, I, II 
and HI turn out to be 2.55, 2.56 and 2.60 respectively. The 
completeness of this selection criterion i.e. the fraction of the 
total number of segments with QPO that would be selected 
turn out to be 0.36, 0.82 and 0.69 for the three states respec- 
tively. The corresponding frequencies of the segments with 
P > Pmax are shown as a function of the normalised count 
rate in Figure 3. As expected, the parallel tracks are repro- 
duced. The histogram of the selected segments for different 
count rates (Figure 4 top panel) shows peaks at the same 
fluxes as in Figure 1. This may indicate that the gaps in 
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Figure 1. The histogram of normalised count rate Rpf for all 
the data sets considered. Four flux states have been identified and 
marked by roman numerals. I: 1225 < Rm < 1625 counts/sec, 
II: 1625 < Rn < 2050 counts/sec. 111: 2050 < Rn < 2550 
counts/sec, IV: Rn > 2550 counts/sec. 
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Figure 2. The histograms of the maximum Leahy normalised 
power Pmax for four flux states. The dotted line is the ex- 
pected probability distribution due to statistical fluctuations 
alone Pns{Pmax)- The best flt parameters, the reduced ^rid 
the number of degrees of freedom, d.o.f, are given. For all proba- 
bility distributions cr = 0.15. 

the frequency-flux plot (Figure 3) maybe due to the uneven 
distribution of fluxes exhibited by this source. However, the 
ratio of the segments with P > Pthre to the total number of 
segments at a flux level (Figure 4 b) also shows peaks at the 
same corresponding fluxes. Thus it appears that there are 
at least three distinct flux levels for this source. The QPO 
phenomenon appears when the system is in one of these flux 
level , while in the intermediate flux levels the phenomenon 
is either weak or absent. 

The histogram of the number of segments with P > 
Pthre as a function of frequency is shown in figure 5, while 
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Figure 3. The variation of the QPO frequency for all segments 
with P > Pthre versus the normalised count rate. 
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Figure 4. The number of segments with P > Pthre> (Top 
panel) and the ratio of Nq to the total number of segments, No 
(Bottom panel) as a function of the normalised count rate (-Rjv)- 



figure 6 shows the average peak amplitude < Aq >, defined 
as Aq = {P — 2)/{0.38N). Altliougli, Aq is a measure of the 
peak of the Lorentzian function, rather than the actual am- 
plitude of the QPO, t he results obtained h e re are similar to 
the re sult obtained bv lMendez' et al. H200lD : lDi Salvo et al. I 
l|20o3), who found that the amplitude is nearly constant till 
/ f» 850 Hz and then falls off. However, as mentioned ear- 
lier these results should be treated with caution since the 
selection criterion makes the QPO sample incomplete. 

For a selection criterion based on flux levels, as has 
been done in this work, the QPO detected may arise from 
different parallel tracks. In principle, each track may have 
a different average amplitude and/or incidence. Such a di- 
chotomy may give rise to a bimodal distribution of QPO 
powers and should be fitted by two probability distribution 
functions. This may be particularly true, when the frequency 
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Figure 5. The number of segments with P > Pthre^ ^Q versus 
the QPO frequency for the three lower flux levels (I, II and III). 

shift between the two tracks is large and the flux range cho- 
sen to define a fiux level is small. However, the flux range 
chosen in this work, is sufficiently large such that the cor- 
responding range of frequencies for each track is large with 
significant overlap. Thus, the distribution of QPO versus fre- 
quency (Figure 5) does not show any bimodal structure for 
states (H) and (HI). For state (I) there may be a double peak 
structure in the distribution. However, for this state the av- 
erage amplitude does not seem to vary much with frequency 
(Figure 6) implying perhaps that the average amplitude is 
nearly same for the two tracks. Fitting the detection dis- 
tribution versus maximum power (Figure 2; top left panel) 
with two probability functions, instead of one, does not lead 
to a reduction in . Instead a degenerate fit where the am- 
plitudes of the two probability distributions are nearly equal 
and the incidence fractions add up to ~ 0.9, is obtained. Any 
difference between the incidence and average amplitude of 
the two parallel tracks may be detectable with the use of 
smaller fiux ranges, but this leads poor statistics due to the 
lower number of detections for each flux level. Alternatively, 
this may be achieved, by a different selection criteria which 
uses some other source property that can clearly distinguish 
and separate out the different parallel tracks. However, it is 
not clear which source property can be used here, since the 
QPO frequency has a weak dependence on hardness ratios 
tPi Salvo et aL II2003I) . 



5 SUMMARY AND DISCUSSION 

Based on the histogram of the number of detections as a 
function of normalised count rate, four distinct flux levels 
have been defined. It is found that the incidence of kHz 
QPO, decreases from near unity at low fiux levels to nearly 
zero at high ones. The variation of frequency with count 
rates reveal roughly three distinct parallel tracks which are 
coincident with the lower three fiux levels. It is however, 
argued that in the intermediate flux levels, the QPO phe- 
nomenon actually disappears or weakens and the parallel 
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Figure 6. The averaged amplitude Aq for segments with with 
P > Pthre^ versus the QPO frequency for the three lower flux 
levels (1,11 and III). 

tracks are not due to uneven sampling of the source at dif- 
ferent intensities. 

This analysis suggests that the QPO phenomenon is 
also related to the long term flux level of the source and in 
particular occurs more frequently as the flux level decreases. 
The existence of an inferred threshold accretion rate, be- 
yond which the QPO phen omenon disappears has been sug- 
gested earlier JCui I |2000D based on short term evolution 
of the source. Here we show that for the long term evo- 
lution there exists such a threshold intensity level. Within 
each long term flux level, it is shown that the QPO weak- 
ens and/or disappears as the count rate increases. Thus the 
pheno menology se ems to be more complex than envisioned 
before. ICui I i)2onnh have proposed that the QPO is produced 
as a result of the interaction of the magnetosphere with the 
Keplerian disk, and at high accretion rates, when the mag- 
netosphere radius is smaller than the last stable orbit, the 
phenomenon disappears. If the magnetosphere radius is the 
radius at which magnetic pressure equals the ram pressure 
of freely accreting material, then this model predicts a sin- 
gle accretion rate threshold above which the phenomenon 
should disappear. Thus this simple model needs to be mod- 
ified to explain the disappearance or weakening of the QPO 
phenomenon for each long term flux level. This may not be 
surprising, because the magnetosphere for this low magnetic 
field systems, probably depends also on the kind of accretion 
flow rather than just on the accretion rate and neutron star 
dipolar field strength. 

The presence of parallel tracks in the QPO frequency 
versus intensity plot, seems to imply that there is more than 
one parameter that determines the QPO frequency. It is 
expected that one variable parameter is the accretion rate 
which should determine the accretion flow and hence QPO 
properties. It is not clear which is the other varying parame- 
ter which gives rise to shifts in the c orrelation leading to the 
parallel tracks. Ivan der Klis I i200l[) showed that this prob- 
lem can be resolved if the QPO frequency depends on both 
the prompt and time-averaged value of a single variable pa- 
rameter. In this scenario, the parallel tracks arise because 



the time-averaged intensity (or accretion rate) is different 
for the different tracks. In the simplest version of this model, 
there is no preferred intensity levels and the distinct tracks 
are obtained due to observational windowing i.e. the tracks 
are not connected due to gaps in observations. However, the 
results obtained in this work suggest that the parallel tracks 
may correspond to definite flux levels and the gaps in the 
plot are not due to under sampling. Thus the above model 
needs to be modified to take into account these observations. 

The QPO phenomenon seems to be correlated to long 
term flux levels on the source in a complex manner, which 
makes simple theoretical interpretations difficult. The prob- 
lem is enhanced because as shown in this work the variation 
of the amplitude of the QPO is affected by the incomplete- 
ness arising from the data selection technique. Further, it has 
been tacitly assumed that in all the flux states, the same 
QPO is being studied. A clearer idea would arise after a 
study of the spectral properties of this system. A systematic 
study of the spectra ( not just the colour which is affected by 
detector response) of the source at different flux levels and 
when the QPO is present or not, would shed light into the 
phenomenon's origin. It may then be possible, to develop 
and test theoretical models regarding the nature of these 
systems. 
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